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Abstract 

A 

An approach is presented for determining the 

minimum weight design of aircraft wing models which h k k 

takes into consideration aerodynamics-structure cou- ’ v ’ 
pling when calculating both zeroth order informa- 
tion needed for analysis and first order information £ 
needed for optimization. When performing sensi- 
tivity analysis, coupling is accounted for by using a ri 
generalized sensitivity formulation. The results pre- f * 
sented show that the aeroelastic effects are calcu- f' 
lated properly and noticeably reduce constraint ap- ' 
proximation errors. However, for the particular exam- s' 

pie selected, the error introduced by ignoring aeroe- 
lastic effects are not sufficient to significantly affect 9l 
the convergence of the optimization process. Trade KS 
studies are reported that consider different structural 
materials, internal spar layouts and panel buckling 
lengths. For the formulation, model and materials p 
used in this study, an advanced aluminum material D ' 
produced the lightest design while satisfying the prob- 
lem constraints. Also, shorter panel buckling lengths 
resulted in lower weights by permitting smaller panel p'» 
thicknesses and generally, unloading the wing skins 
and loading the spar caps. Finally, straight spars re- i 
quired slightly lower wing weights than angled spars Po 


Nomenclature 

Scalar quantities are typed in lower case and 
normal script, column vectors are lower case and bold 
scnpt, matrices are upper case and bold script. Also 
the superscript l refers to one of the flight conditions 
while t indicates the transposed of a matrix or vector 
Finally, four different grids are used for reference of 
vector quantities: 1) the aerodynamic pressures and 
forces are calculated on a 15x15 grid overlaid over 
the aircraft planform, ii) the concentrated masses and 
inertia forces are found on an uneven inertia grid, iii) 
both pressure and inertia grids are combined on a 
total gnd, iv) constraints are calculated on 5x5 grids 
overlaid over the trapezoidal panels making up the 
aircraft configuration. 


P' 


n l ,v l ,w l 


u 

u', v', w' 


W' 


W e 


w. 


diagonal matrix containing 
reference areas on pressure grid 

elastic displacement modes in x,y 
and z directions 

aerodynamic pressure loads, on 
pressure grid 

inertia loads, on inertia grid 

total point loads, on total grid 

vector of constraints 

constraint function 

Kresselmeier-Steinhauser 
envelope constraint 

total pressures, on pressure grid 

pressures due to camber, on 
pressure grid 

pressures due to rigid camber, on 
pressure grid 

pressures due to a unit 
angle-of-attack, on pressure grid 

pressures due to elastic 
displacement modes, on pressure 
grid 

vertical load factor 

polynomial approximation of elastic 
displacements 

vector of 1 

generalized displacements 
aircraft empty weight 
vector of lumped empty weights 
fuel weight 
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vector of lumped fuel weights 

aircraft gross weight 

chordwise coordinate 

center of gravity longitudinal 
position, aircraft with fuel 

center of gravity longitudinal 
position, aircraft without fuel 

fuel center of gravity longitudinal 
position 

center of pressure 

vector of longitudinal positions of 
concentrated load application 
points 

independent structural design 
variables 

spanwise coordinate 

dependent (behavior) variables 
calculated by aerodynamics and 
used in structures 

dependent (behavior) variables 
calculated by structures and used 
in aerodynamics 

dependent (behavior) variables 
calculated by structures and not 
used anywhere in another 
subproblem 

trim angle-of-attack 

user specified coefficient in the KS 


Introduction 

This paper reports on a study to formally carry 
out the integration of aerodynamic and structural 
analyses in the minimum weight design of a super- 
sonic transport aircraft wing. This is a necessary 
first step in carrying out the complete integration of 
aerodynamic, performance and structural analyses 
in the optimization of the shape and structural de- 
sign of a wing for optimum performance. The gen- 
eral approach to solve the latter problem has been 
discussed by Barthelemy et a!> and will not be re- 
peated here. Instead, the focus is on the details of 
the formulation of the static aeroelastic problem, on 
the implementation of the optimization procedure and 
also, on the latest technical results. 

The capability to perform minimum weight struc- 
tural design under static aeroelastic constraints has 


existed for a long time. In most applications, how- 
ever, the aeroelastic coupling is only accounted for 
while performing system analysis, that is, while com- 
puting zero-th order information like loads, displace- 
ments, stresses and strains. Specifically, the cou- 
pling effect is ignored while performing system sensi- 
tivity analysis, that is, calculating first-order gradient 
information necessary for optimization. In this paper, 
the coupling is explicitly accounted for in both analy- 
sis and sensitivity analysis. The tool used to explic- 
itly account for the coupling between the disciplines 
is the Generalized Sensitivity Equations discussed at 
length by Sobieski 2 . 

This study is part of HiSAIR (High-Speed Aircraft 
Integration Research), a research effort at the NASA 
Langley Research Center. As discussed by Dollyhigh 
and Sobieski 3 HiSAIR is an ”... organized program to 
strengthen multidisciplinary aircraft analysis, design 
and optimization ...”. It “... cuts across organizational 
discipline lines of the Center and uses the High- 
Speed Civil Transport as a research focus.” Coen 
et al. A give a thorough activity-by-activity review of 
the status of HiSAIR. The current vehicle focus for 
HiSAIR is a civil transport aircraft designed to carry 
252 passengers for 6500nm at a Mach number of 2.4 
at a mid-cruise altitude of 63000ft. 

The paper begins with a discussion of the prob- 
lem formulation, reviewing the specifics of the struc- 
tural, aerodynamic, and aeroelastic analyses and 
sensitivity analyses. Then the optimization procedure 
is discussed. Finally a section describes current re- 
sults; it includes a discussion of the trade studies that 
consider different materials and internal spar layouts. 

Problem Formulation 


Structural Analysis and Modelling 



alent laminated plate formulation implemented in 
program ELAPS (Giles 56 ), a code based on the 
Rayleigh-Ritz method. It can model wing skin as well 
as rib and spar caps. In ELAPS, the analysis process 
begins with a wing described by planform geome- 
try, thickness, camber, and material information and 
modelled by several trapezoidal plates. ELAPS then 
proceeds by assuming that the three displacement 
fields for the wing can be represented by polynomi- 
als in the chordwise (z) and spanwise (y) coordinates. 
For each flight condition /: 

u'(z,y) = b£(z,y)u' 

v'( x >y) = b v( a; .i/) v( W 

w'(x,y) =b^(z,y)w' 

where the entries of vectors b u ,b v ,b w are referred 
to as elastic displacement modes; they are products 
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of the type x n y m . The coefficients u',v',w' are 
generalized displacements. 

Figures 1 and 2 depict the aircraft models con- 
sidered in this study; they differ by whether the spars 
are straight and swept or cranked so as to be orthog- 
onal to the fuselage axis. Each model is made of 1 1 
separate plates, one of which (plate 11) constitutes a 
model of the fuselage scaled to match known fuse- 
lage weight data. The model has 4 spars and 10 ribs 
whose caps can be modelled by ELAPS. The bound- 
ary conditions specify symmetry of the displacements 
about the centerline. In addition, the model is sup- 
ported in pitch by two springs; however, as explained 
later, the aerodynamic and inertia loads are balanced 
in each flight condition, so that the reactions in the 
springs are zero or nearly so. The wing skins are 
sandwich panels with ,75in thick honeycomb core. 



Figure 1 Swept spar model 
definition, plate and spar numbering 



Figure 2 Cranked spar model 
definition, plate and spar numbering 

Earlier comparisons between finite element and 
ELAPS results for similar models have pointed to 
the difficulty to capture stresses in models with com- 
plex planform, thickness and camber distributions 
adequately. While it would be desirable in such 
cases to use higher-order polynomials, the current 


implementation of ELAPS, based on simple polyno- 
mial displacements, does not permit it. Therefore, 
to obtain reasonable stress distributions, the current 
model neglects the camber of the wing and calculates 
stresses, assuming a symmetric airfoil. The highest 
power appearing in the elastic displacement modes 
are x 4 y 5 for u' and v', and x A y 6 for w'. This results 
in 25 generalized degrees-of-freedom (dofs) for dis- 
placements u' and v' and 30 dofs for displacement 
w'. The whole model therefore has 80 dofs. 

Given a selection of b’s, and point forces f/, 
ELAPS calculates the generalized displacements. 
From those displacements, it finds stresses and 
strains in the wing covers. These are evaluated on 
upper and lower wing skins, on a 5x5 grid on each 
individual plate. Conventional constraints are then 
evaluated for the three plane strains, the von Mises 
stress as well as orthotropic panel buckling. The 
buckling constraint is an orthotropic buckling equa- 
tion given for a simply supported plate by Sobieski 7 ; 
the critical shear load of an infinite length strip is taken 
from the DOD/NASA Advanced Composite Design 
Guide 8 . 

ELAPS also finds strains and stresses evaluated 
at five equally spaced stations on each cap. Stress 
and strain constraints are calculated at those stations. 
Finally, ELAPS also evaluates the structural weight. 

Only plates 1-10 are retained when calculating 
constraints as it is felt that the fuselage stresses are 
not representative of actual stresses. This leaves 500 
each of the plate constraints and 460 each of the 
cap constraints. This makes for a very large num- 
ber of constraints for the optimizer to handle. To re- 
duce that number, an envelope function developed 
by Kresselmeier and Steinhauser 9 is used to com- 
bine together all constraints of the same type for a 
given flight condition; this leaves 5 plate envelope 
constraints and 2 cap envelope constraints for each 
flight condition. For a set of individual constraint func- 
tions gj, the envelope constraint is defined by: 

K S = 1 tn exp (pgj ) j (2) 

where p is user-selected. Parameter p must be care- 
fully selected (see Barthelemy and Riley 10 , for exam- 
ple) the larger the value of p, the closer the envelope 
is to the initial constraints boundaries. On the other 
hand, excessively large values of p cause the enve- 
lope to have nearly discontinuous gradients. 

The wing panel face sheet thicknesses at each 
panel comer are taken as design variables and the 
skin thickness is assumed to vary linearly between 
the comers. The spar cap cross-sectional areas at 
the plate boundaries are design variables and the 
cross sectional areas vary linearly between plate 
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boundaries. Finally, the rib cap cross-sectional ar- 
eas are variables and each rib cross-sectional area 
is constant between leading edge and trailing edge 
spars. The upper and lower wing covers are identi- 
cal. There is a total of 44 design variables describing 
the model, they are contained in vector x,. 

Aerodynamic Analysis and Modelling 

Aerodynamic analysis is performed with the 
WINGDES computer program from Carlson and 
Walkey 11 . WINGDES computes subsonic and super- 
sonic pressure distributions on a 15x15 grid overlaid 
on the aircraft planform. In both modes, it uses a lift- 
ing surface formulation. For each flight condition, the 
linearity of pressure distributions enables the defini- 
tion of pressure due to angle-of-attack (zero camber) 
p' tt and pressure due to camber (zero angle-of-attack) 
pj.. The pressure due to camber is due to the rigid 
camber (and twist) built into the wing as well as to 
the elastic deformations resulting from the airloads. 
The pressure due to rigid camber is defined as p„, 
while the pressure due the transverse displacement 
«''(*. y) = ,y) w' can be found by superposition 

as 

p' e = P ,7 V (3) 

where P' is a matrix, each column of which is a 
vector of pressures due to one of the transverse 
displacement modes. It must be noted that vectors 
Po>Po and matrix P' are fixed for any given flight 
condition and wing planform. 


Integrated Static Aeroelastic Analysis 

The object of integrated aeroelastic analysis is to 
find the airloads and resulting elastic wing displace- 
ments which result in a trimmed aircraft in the various 
flight conditions. Those conditions are summarized in 
Table 1. Trimming is effected by selecting the angle- 
of-attack of flight that balances the total vertical load 
on the aircraft as well as the fuel distribution that bal- 
ances the pitching moment. Since the aircraft is tail- 
less, the later condition requires that the aircraft cen- 
ter of pressure and center of gravity are collocated. 

The given input geometry corresponds to the 
cruise condition (flight condition superscript / = 1) . 
The jig shape is obtained by removing cruise elastic 
displacements w‘ from the input camber. In any off- 
cruise flight condition, the elastic deformations from 
cruise shape are given by w* -w 1 . 

The total pressure due to camber can now be 
written: 

p' =p' +P' T (w'-w‘) (4) 

so that the total aerodynamic pressure, including 

elastic effects is given by 

p' = P c + <*Va (5) 

and, the total aerodynamic forces are then: 

f a = Ap' (6) 

where A is a diagonal matrix whose entries are the 
areas associated with the reference pressure points. 


Table 1 Flight conditions description 


Load case 

Load factor 

Mach 

number 

Altitude (ft) 

Dynamic 

pressure 

(lb/ft 2 ) 

Fuel weight (lbs) 

1, Mid-cruise 

1.0 

2.4 

63175 

535 

85544 

2, Transonic climb 

1.0 

1.2 

29670 

644 

167016 

3, Reserve cruise 

1.0 

0.9 

42480 

199 

18446 

4, High-Speed pull-up 

2.5 

2.4 

56949 

723 

149830 

5, Low-speed pull-up 

2.5 

0.6 

10000 

367 

174185 
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The vertical force equilibrium equation is used 
to find the trim angle-of-attack. Neglecting rotary 
inertias: 

n‘w' s = u T f' = u T A(p' + a'pj.) (7) 

where all the entries of vector u are 1. 


The aircraft is trimmed in pitch by moving fuel 
among the different tanks so that the center of gravity 
coincides with the center of pressure x[ g = x' cp . The 
center of pressure location is given by: 

f/T 

l 1/1 An 


X c P = 


f' T u 


The total center of gravity location is given by: 


x cg w g = x r 


.U>e + X c 


( 8 ) 


(9) 


Pitch equilibrium determines a target position for the 
fuel center of gravity that is given by 


x cgJ - 


X CS W 9 ~ X 


cge 1 


(10) 


To obtain that target position, the existing fuel is dis- 
tributed in the tanks, starting from those closest to 
the target position, so that their moment with respect 
to that position remains exactly zero. With the cur- 
rent tank layout, this process has proven generally 
feasible and extremely fast. It must be noted that 
while this approach is generally correct for straight 
and level flight conditions, it is quite approximate for 
pull-ups since i) rotary inertias are neglected and ii) 
it would not be possible to move fuel fast enough to 
trim the aircraft in such a condition. 


The inertia loads are given on the inertia grid, 
they become 

f' = n'(w e +w') (11) 


And the total vector of concentrated forces be- 
comes, by concatenation, 

f t ' = (f ‘ T ,Cf (12) 


Integrated Aeroelastic Sensitivity Analysis 

The formulation outlined above shows coupling 
of aerodynamic and structural disciplines. Integrated 
sensitivity analysis can be applied to obtain total sen- 
sitivity derivatives that account for the coupling. This 
requires clearly identifying the disciplines interacting 
in the coupled system and delineating the boundaries 
between those disciplines. 

The coupled set of equations for the aeroelastic 
system described above can be written functionally 
as follows: 


y as — Yax (y sa ) 

y»a = y a a(Xg,y M ) ( 18 ) 

y*a = y« a (x.,yas) 


Equation 13 identifies the relationships between the 
outputs of a discipline and its inputs. These inputs 
may be either independent variables or dependent 
variables calculated by another discipline. In this no- 
tation, the first subscript identifies the discipline where 
the dependent variable is calculated, the second sub- 
script, where it is used. When the subscript is re- 
peated, the dependent variable is calculated but not 
used in any other discipline analysis; it is used in the 
design process. Figure 3 gives a graphic description 
of the problem, identifying the disciplines as well as 
the independent and dependent variables. 
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Figure 3 Problem functional description, 
independent and dependent variables definition 


For given flight conditions and aircraft planform 
and camber distributions, the analysis begins with 
generation of vectors p' a , pj, and matrix P'. Then, for 
each flight condition, an iterative loop ensues that first 
finds the pressure due to camber, then the trimmed 
angle of attack, fuel distribution as well as the total 
pressure distribution and the total loads on the wing. 
Structural analysis is performed to calculate a new 
elastic camber and the process is restarted. The 
iterative loop continues until convergence. 

This procedure was validated by comparing its 
predictions with those from the TSO program 12 , for 
a simple trapezoidal swept wing trimmed in angle of 
attack only. In a 2.5g pull-up at a Mach number of 
.45, the approaches yielded differences of 5% in rigid 
and elastic trim angles of attack and 3% in wing tip 
displacement. 


Vector y M contains the dependent variables cal- 
culated in the aerodynamic discipline and needed in 
the structural discipline for calculation of displace- 
ments, strains and stresses. In principle, this would 
include vectors of aerodynamic pressures in the dif- 
ferent flight conditions. However, because of the lin- 
ear representation retained for the wing pressure, and 
of the fact that vectors p‘ a , p‘ 0 and matrix P' remain 
constant throughout the design, this information is 
limited to the amplitudes of the pressure modes in 
Eq. 4, that is the vectors of generalized displace- 
ments measured with respect to cruise shape. 

y M = 

f(w a — w 1 ) T ,(w 3 — w l ) T ,(w 4 — w 1 ) T , (w 5 — w 1 ) 7 ^ 

(14) 
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Vector y sa contains the dependent variables calcu- 
lated in the structural discipline and needed in the 
aerodynamic discipline for calculation of pressures. 
This includes the generalized displacements in the 
different flight conditions. 

y aa = (w lT ,w aT ,w 3T ,w 4T ,w 5T ) T (15) 

Finally, y sa , the vector of dependent variables calcu- 
lated in the structures discipline and not needed by 
other disciplines contains the aircraft gross weight at 
cruise and the envelope constraints in each of the 
load cases. 

y«. = (wg, g lT .g aT ,g 3T ,g 4T .g 5T ) T (16) 


The generalized sensitivity equations read: 


I 

-fc °1 


/ 0 

fly.» 
oyl 7 

7 

0 

(fr = 

[p 

Sy*. 

0 

7. 

Vte/ 

' dx m 


(17) 


block in the decomposition. Rather, each block can 
be thought of as a computational module performing 
a set of related calculations. Each classical disci- 
pline can then comprise several blocks. It turns out 
that there can be several valid decompositions of the 
problem. The elements to take into account when 
identifying the blocks include i) what are the natural 
computational boundaries (it is impractical to define 
a boundary in the middle of a computer program), ii) 
what are the modules for which partial derivatives are 
already available or easy to obtain, iii) what are the 
dependent variables for which total sensitivity infor- 
mation is required. 

For example in this problem, one could have 
added a third discipline dealing with the balance of 
the aircraft. This would have added to the complexity 
of the decomposition but would have generated more 
information, including the sensitivity of the aircraft 
trim angle of attack and center of gravity position to 
changes in the design variables. 


Optimization 


This equation gives the sensitivity derivatives of the 
dependent variables accounting for the coupling be- 
tween the disciplines: the total derivatives (d(.)/d(.)) 
as a function of the sensitivity derivatives ignoring the 
coupling between the disciplines: the partial deriva- 
tives (d(.)/d(.)). 

For example, dy^/dy^ is the matrix of deriva- 
tives of the amplitudes of the pressure modes with 
respect to the amplitudes of the displacement modes; 
using Eqs 14 and 15, it is quite simple to calculate 
and is given by 


<9ya» 

dy.a 


■-/ I 0 0 o- 

-I 0 7 0 0 
-7 0 0 /0 
.-7 0 0 0 7. 


(18) 


Also, <9y. a /dyas and dy M /dy M are matrices of 
derivatives of the elastic deformations, and the 
weight and stresses, respectively, with respect to 
the amplitudes of the pressure modes. Finally, 
dy«a/dx s and dy M /dx s are matrices of derivatives of 
the elastic deformations, and the weight and stresses, 
respectively, with respect to the structural design vari- 
ables. The last four derivative matrices are calculated 
by forward finite difference. In this application, this 
detailed formulation permits one to account for elas- 
tic load redistribution while calculating the gradients 
of the dependent variables, including the constraints. 

There is a certain degree of arbitrariness in iden- 
tifying the disciplines involved and in deciding where 
the boundary lines have to be drawn. In fact, there 
is no specific requirement that the calculations be- 
longing to a given discipline be confined to a single 


Figure 4 presents a graphic description of the 
generic optimization capability developed in this 
study. It is a SUN engineering workstation based 
system currently implemented to handle 5 disciplines 
with up to 100 independent variables and 500 depen- 
dent variables. The system is designed to provide for 
user intervention at any point in the design process. 
Alternately, the system can operate completely auto- 
matically. It proceeds in design cycles, each requiring 
full analysis and sensitivity analysis of the problem. 
Within each cycle, it is possible to change such things 
as the type of problem approximation, the type of op- 
timization algorithm used, the combination of depen- 
dent and independent variables optimized, the move 
limits for approximations. Each change of one or sev- 
eral of these parameters produces a different design 
alternative. 
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Figure 4 Integrated design system 
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The system is built around the commercial pack- 
age OPTDES 13 which offers several optimization al- 
gorithms. Those used in this study are linear pro- 
gramming, sequential linear programming, method of 
centers, generalized reduced gradients, and sequen- 
tial quadratic programming. Since analyses and sen- 
sitivity analyses are quite expensive, OPTDES op- 
timizes a sequence of approximations to the actual 
design problem. These approximations are all based 
on zeroth and first order information of the depen- 
dent variables and include linear, reciprocal, convex 
and the two-point approximation of Fadel et a/. 14 . For 
the results generated in this paper, the combination 
of two-point approximation and sequential linear pro- 
gramming proved the best. 

Because of the use of approximations, move 
limits must be imposed to guarantee convergence of 
the optimization process. However, move limits may 
prevent design progress from an infeasible starting 
point. To avoid this problem, a constraint relaxation 
procedure similar to that discussed by Barthelemy 
and Riley 10 is implemented. 

To provide an audit trail for the design process 
and allow for restart from any design cycle, critical 
optimization information is stored primarily in the Re- 
lational Information Manager (RIM 15 ), a commercial 
relational database management system. Cycle in- 
formation retained includes initial values and upper 
and lower bounds on the independent and depen- 
dent variables. Because of its potential size, cycle 
gradient information is kept in conventional file for- 
mat. Design alternative information retained includes 
final independent and dependent variables for each 
alternative design within each cycle. 

Each design cycle begins with system analysis 
and subsystem sensitivity analysis. This step can be 
conducted with any existing analysis package on any 
computer or distributed system of computers. Each 
discipline produces one file containing its own anal- 
ysis and sensitivity analysis information. This infor- 


mation is then processed to produce the generalized 
sensitivity equations (GSE, Eq. 17) and relevant data 
is stored in the RIM database and the gradient files. 
Once optimization is completed, the user may inter- 
actively query the database and track graphically or 
in tabular output any combination of independent or 
dependent variables. The user may also gauge the 
accuracy of the approximations selected by compar- 
ing analysis results predicted with those obtained af- 
ter reanalysis, and then decide to produce more de- 
sign alternatives within the current cycle or to initiate 
a new cycle using as the starting design any of the 
design alternatives generated previously. 

Results 

Optimization Results 

The procedure described in this paper was used 
to find the minimum weight design of models with 
two different spar layouts, two different materials and 
three different panel buckling lengths. In changing 
the buckling length, however, the weight of the spar 
and rib webs was assumed to remain constant. Ta- 
bles 2, 3, 4 and 5 give summary results for optimiza- 
tion. In these tables, the structural weight includes 
wing skins and rib and spar caps for both covers of 
one wing. The active variables are those not locked 
against their lower bounds; lower bounds are set at 
.004 ft for the skin thicknesses, .005 ft 2 for the cap 
cross-sectional areas. For the plates, the constraints 
are identified as v' for von Mises, b' for buckling, x' 
for chordwise strain, y' for spanwise strain, xy' for 
shear strain. For the caps, e' is used for strain, s' 
for stress. Also, the superscript / stands for the flight 
condition. The constraints are written so as to range 
between -1. and .0 when satisfied. In the tables, the 
constraints are identified as active if between -.02 
and .0 and as nearly active (as indicated by paren- 
theses) if between -.1 and -.02. 


Table 2 Summary of final designs, all-titanium wing, swept spars. 


Buckling Length (ft) 

3.0 

2.5 

2.0 

Structural Weight (lbs) 

33058 

26525 

24698 

Active Variable Numbers 

plates 

3, 4, 5, 6, 7, 8 

6, 7 

6, 7 

spars 


25, 29 

24, 25, 28, 29 

| Active Constraints (superscript indicates load case) 

plates 

b 5 

v 2 , b 5 

v 2 , v 5 , (xy 2 ) 

spars 



(s 2 ), (s 5 ) 
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The active thickness variables are the thick- 
nesses at the comers of the plates corresponding to 
the wing box (Figs. 1 and 2), that is: plates 2 (vari- 
able number 3, 4, 5, and 6), 3 (5, 6, 7, and 8) and 
4 (7, 8, and 9). For the swept spars, the active vari- 


ables are for spars 2 (24, 28, and 32) and 3 (21, 25, 
and 29); for the cranked spar case, spar 1 variables 
(19 and 23) become active also. Rib cap variables 
never stay active at the optimum. 


Table 3 Summary of final designs, all-aluminum wing, swept spars 


Buckling Length (ft) 

3.0 

2.5 

Structural Weight (lbs) 

25791 

23295 

Active Variable Numbers 

plates 

4, 6, 7, 8, 9 

6, 7.8 

spars 

25, 29 

21, 24, 25, 28, 29, 32 

Active Constraints (superscript indicates load case) 

plates 

v 2 , b 5 

(v 2 ). (xy 5 ) 

spars 


(s 2 ) 


Table 4 Summary of final designs, all-titanium, cranked spars 


Buckling Length (ft) 

3.0 

2.5 

2.0 

Structural Weight (lbs) 

34931 

28947 

25130 

Active Variable Numbers 

plates 

5, 6, 7, 8, 9 

5. 6, 7, 8, 9 

5, 6,7 

spars 

26 

25,28 

19, 23, 24, 28 

| Active Constraints (superscript indicates load case) 

plates 

v 2 , b 5 , (xy 2 ) 

v 2 , b 5 , (xy 2 ) 

v 2 . v 5 , (y 5 ), (xy 2 ) 

spars 



s 5 , e 5 


Table 5 Summary of final designs, all-aluminum, cranked spars 


Buckling Length (ft) 

3.0 

2.5 

2.0 

Structural Weight (lbs) 

28391 

24235 

22779 

| Active Variable Numbers 

plates 

5, 6, 7, 8 

5. 6, 7, 8 

5, 6, 7, 8 

spars 

25, 28,. 29 

19,25, 28, 29 

25, 28, 29 

| Active Constraints (superscript indicates load case) 

plates 

v 2 , b 5 , (xy 2 ) 

v 2 , b 5 , (xy 2 ), s 5 

v 2 , (y 5 ), (xy 2 ) 

spars 



(s 2 ), s 5 , (e 5 ) 


8 
































For each material and spar layout combination, 
the same trend develops with panel buckling length. 
The optimum design with the largest panel length 
is designed primarily by the buclding constraint in 
the low-speed pull-up; notably no cap constraint is 
active and few cap variables are active . As the 
panel buckling length is decreased, buckling criticality 
decreases, thinner skins can be used and the loads 
are transferred to the spars. In the process, more 
cap variables and constraints become critical. This 
trend is a consequence of the facts that i) spars 
and skins use the same material, and ii) because 
of the geometrical layout of the wing sections and of 
the definition of skin thickness variables, the spars 
offer more bending resistance for a given amount of 
structural weight. Therefore, if panel buckling is not 
an issue, it is more efficient to put the material in the 
spar caps. 

Figures 5, 6 and 7 demonstrate the trend graphi- 
cally for the all-titanium design with swept spars. Fig- 
ure 5, traces the buckling constraint in flight condition 
5, the low-speed symmetric pull-up. At the starting 
point, the constraint is violated on both the inboard 
and outboard sections. The 3ft buckling length op- 
timum design has the buckling constraint active on 
both the inboard and the outboard wing sections. As 
the buckling length is reduced, the constraint remains 
only active on the inboard section of the wing and 
then becomes inactive for a 2ft buckling length. In 
contrast, Fig. 6 traces the von Mises criterion in flight 
condition 2, the transonic climb. Initially violated in 
the outboard section of the wing, the constraint is 
satisfied at the largest buckling length, and becomes 
active or nearly active for both inboard and outboard 
panels at the smallest buckling length. 

Figure 7 shows the various optimum thickness 
distributions. Each has a maximum at the break be- 
tween the inboard and outboard wing sections, over 
the trailing edge spar. The maximum thickness de- 
creases noticeably with the buckling length. While 
none of the active constraints are ever located at the 
wing break, the fact that the maximum thickness oc- 
curs at that location is directly related to the definition 
of skin thickness variables and skin thickness vari- 
ation which guarantees a maximum thickness at a 
plate corner. Since the active constraints peak either 
in the inboard wing section or the outboard wing sec- 
tion, it is logical that the maximum thickness occurs 
at the break section, which defines the boundary be- 
tween plates 2 and 3. 

Additional comparison of Tables 2, 3, 4 and 5 
shows the swept spar designs to be very slightly 
lighter than the cranked spar ones. The advanced 
aluminum designs are lighter than the titanium ones, 
more so for the buckling dominated case; this is a 
direct consequence of the better specific stiffness of 
aluminum. 



B) PANEL BUCKLING LENGTH 3 FT 



C) PANEL BUCKLING LENGTH 2.5FT 


BUCKLING 



D) PANEL BUCKLING LENGTH 2 FT 


Figure 5 Effect of panel buckling length 
on buckling criterion, upper skin, flight 
condition 5, all-titanium swept spar model 

Optimization for given material, spar layout and 
buckling length proceeded with between 20 and 50 
cycles per design, each cycle lasting 1 hour and 40 
minutes on a combination of one IRIS and two SUN 
engineering workstations. The longest runs were re- 
quired with the largest buckling length where buckling 
was an active constraint. As shown on Fig. 5. two 
different plates have critical or near-critical individual 
buckling constraints in that case. The envelope con- 
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A) INITIAL DESIGN 



B) PANEL BUCKLING LENGTH 3 FT 



C) PANEL BUCKLING LENGTH 2.5FT 


VON MISES 



D) PANEL BUCKLING LENGTH 2 FT 


Figure 6 Effect of panel buckling length 
on von Mises criterion, upper skin, flight 
condition 2, all-titanium swept spar model 

straint combining the individual buckling constraints 
in a single constraint for optimization (Eq. 2) exhibits 
a quickly-varying, nearly discontinuous gradient and 
oscillates between two designs, each dominated by 
one of the individual buckling constraints. To allevi- 
ate the problem, optimization was started with a low 
value of p, which was eventually increased to obtain 
sufficient accuracy. In the problems described here, 
p ranged from an initial value of 10 to a final value 


A) PANEL BUCKUNG LENGTH 3 FT 



B) PANEL BUCKLING LENGTH 2.5 FT 


THICKNESS 



C) PANEL BUCKUNG LENGTH 2 FT 

Figure 7 Effect of panel buckling 
length on optimum skin thickness 
distribution, all-titanium swept spar model 

of 100. At the same time move limits starting at 20% 
and ending at 5% were used to maintain the quality 
of the approximations. 

Sensitivity Analysis and 
Approximations Validation 

It remains to demonstrate the validity and use- 
fulness of including aeroelastic effects when calcu- 
lating gradients. Figure 8 compares linear approxi- 
mations with and without aeroelastic effects included 
in the gradient calculations for the cap stress con- 
straint in the high-speed pull-up. The reanalyses 
were performed for 6 designs located on a straight 
line between the cycle 1 and cycle 3 designs of the 
all-titanium model with swept spars and 3ft buckling 
length. It is shown that including the aeroelastic ef- 
fects noticeably improves the quality of the approxi- 
mations, reducing the error in the slope of the approx- 
imation by about a factor of 4. Figure 9 compares the 



quality of the various approximations available in this 
study for the envelope buckling constraint in the same 
flight condition. It shows the linear approximation to 
be of lesser quality, while all three other approxima- 
tions are comparable. 



Figure 8 Cap stress envelope constraint 
for high-speed pull-up, all titanium swept 
spar design, 3ft panel buckling length 



for high-speed pull-up, all titanium swept 
spar design, 3ft panel buckling length 

The effect of accounting for aeroelasticity on the 
rate of convergence of optimization is shown on Fig. 
10 where the all-titanium, swept spar design is re- 
run. As can be seen, there is no clear benefit on 
convergence of including these effects. Both opti- 
mizations converge to essentially the same design, 
the final weights being within 300lbs. This is con- 
sistent with the fact that studies show that this model 



DESIGN CYCLES 

Figure 10 Optimization convergence, all 
titanium swept spar design, 2ft panel length 

does not exhibit radical aeroelastic load redistribution 
effects, particularly in flight conditions 2 and 5, which 
are critical in this design exercise. So, while including 
the aeroelastic effects improves the quality of the gra- 
dients and of the resulting approximations, the error 
introduced by neglecting these effects is not sufficient 
to alter the overall performance of the optimization 
algorithm. For this design, this indicates that some 
simplifications can be made in the formulation and 
that the dy as /dy. a , dy sa /dy at) dy„/dy a , , dy, a /dx, 
derivatives (Eq. 17) need not be calculated. It must 
be noted, however, that this conclusion is problem- 
dependent and may not be applicable to other design 
problems, other level of analysis detail or even other 
flight conditions. 


Concluding Remarks 


The results discussed in the paper show that 
the aeroelastic effects on sensitivity information are 
calculated properly by the generalized sensitivity ap- 
proach. In addition, accounting for these effects re- 
duces noticeably approximation errors. However, for 
the particular example selected, the errors introduced 
by ignoring the aeroelastic effects are not sufficient to 
significantly affect the convergence of the optimiza- 
tion or the final solution. 

Although the optimization procedure generally 
performs satisfactorily, oscillations occasionally ap- 
pear when two of the same individual constraints are 
nearly active simultaneously for the same flight con- 
dition. In this case, envelope constraint gradients are 
nearly discontinuous and convergence is noticeably 
slowed. This problem is alleviated by adjusting the 
envelope function. At the beginning of optimization, 
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the envelope is kept further away from the bound- 
aries, thereby maintaining smooth gradients; as the 
optimum is neared, the envelope is brought closer to 
the original constraint boundaries. 

Results are given that consider different struc- 
tural materials, internal spar layouts and panel buck- 
ling lengths. For the formulation, model, and mate- 
rials used in this study, an advanced aluminum ma- 
terial generated the lightest design while satisfying 
the problem constraints. Also, lower panel buckling 
lengths resulted in tower weights by permitting tower 
panel thicknesses and generally, unloading the wing 
skins and loading up the spar caps. Finally, swept 
spars generated slightly lighter designs than cranked 
spars. 
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Appendix A Material Properties, 
Allowables, Safety Factors 


Table A1 gives the material properties and al- 
lowables used in this study. Fatigue allowables were 
used in the first three flight conditions. Table A2 
shows the safety factors selected in this problem. In 
general, a factor of 1.5 is used in the high-g toad 
conditions. For the one-g flight conditions, where al- 
lowables are knocked down significantly to prevent 
fatigue of the material, a factor of 1.0 is used. For 
the buckling conditions, an additional safety factor of 
1.5 is used. 
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Table A1 Material properties and allowables 


Material 

Titanium 

Advanced Aluminum 

E (psi) 

16. 10 6 

12. 10 6 

G (psi) 

6. 10 6 

4.78 10 6 

V 

.332 

.318 

p (pci) 

.163 

.105 

cr x = <T y , yield (psi) 

±91.8 10 3 

±61.5 10 3 

&xy , yield (psi) 

±34.5 10 3 

±23.3 10 3 

e* = fy , yield 

±5.74 10‘ 3 

±5.13 10' 3 

fry, yield 

±5.74 10- 3 

±4.82 10' 3 

a z - cr y , fatigue (psi) 

±25.0 10 3 

±15.0 10 3 

c T xy , fatigue (psi) 

±9.38 10 3 

±5.69 10 3 

e x = c y , fatigue 

±1.56 10 3 

±1.25 10- 3 

fry , fatigue 

±1.56 10‘ 3 

±1.19 10' 3 


Table A2 Safety factors used in the various constraints 


Load Cases 

1,2,3 4,5 

Stress, Strain 
Buckling 

1.00 1.50 

1 .50 2.25 
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